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ABSTRACT 

The Anglo- Australian Planet Search has now accumulated 12 years of radial- 
velocity data with long-term instrumental precision better than 3 ms _1 . In this 
paper, we expand on earlier simulation work, to probe the frequency of near- 
circular, long-period gas-giant planets residing at orbital distances of 3-6 AU - 
the so-called "Jupiter analogs." We present the first comprehensive analysis of 
the frequency of these objects based on radial-velocity data. We find that 3.3% of 
stars in our sample host Jupiter analogs; detailed, star-by-star simulations show 
that no more than 37% of stars host a giant planet between 3-6 AU. 

Subject headings: planetary systems - techniques: radial velocities 



1. Introduction 

Recent discoveries of ever lower-mass planets have received a great deal of publicity. 
Equally important, however, are discoveries of long-period planets. A long-standing question 
in astrophysics is "How common are planetary systems like our own Solar system?" The 
so-called "Jupiter analogs," with orbital periods P > 10 years and velocity amplitudes K ~ 
10 ms _1 , represent another means for probing the frequency of systems with architectures 
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similar to our Solar syste Such planets are now within the reach of the longest-running 
radial-velocity planet-search programs. A critical part of addressing the key question of 
Solar system frequency is assessing the selection effects at work in the regime of long-period, 
Jovian planets. 

Microlensing observations are beginning to be effective in detecting and constraining the 
population of Jupiter-like planets. The gravitational microlensing technique is well-suited 
to detecting planets which have not experienced significant inward migration, typically at 
separation s a > 3 AU (cf. Fig . 9 from Sumi et al. 2010). The discovery of a Jupiter/Saturn 
analog by iGaudi et al.l ( 120081 ) illustrates the power of the microlensing technique to detect 
planets which are curre ntly beyond t he re ach (in both observation time and precision) of 
Doppler measurements. iGould et al.l (120 10f ) presented the first estimate of the frequency of 
planets beyond the snow line, based on a sample of 6 microlensing planet discoveries. 

Previously published analyses of radial- velocity data sets have been able to place upper 



l imits on substellar and planetary companions from radial-velocity surveys. iMurdoch et al. 



( 119931 ) determined detection limits for the Mt. John radial- velocity program by adding the 
program's mean velocity error (65 m s -1 ) to the signals of planets in circular orbits. Planetary 
signals were then recovered by the periodogram and F-test methods, and those planets for 
which 95% of phases were recovered with false-alarm probability (FAP)<1% were considered 
detectable. They were able to exclude planets with m sin % >10 Mj up with periods less 
than 2000 days, and b rown dwarf companions (10-40 Mj up ) with P < 8.2 years. Similarly, 
Cumming et al.l ( 119991 ) computed detection limits from the Lick planet search data by 1) 
noting the highest peak z max in the periodogram for each target, and 2) generating simulated 
data sets with sinusoidal (circular-orbit) signals and finding the velocity amplitude K for 
which 99% of signals had power exceeding z max . They achieved a detection lim it of 2 
_1 for companions with a<5 AU. The 12-year CFHT survey of I Walker et al.l ( 119951 ). 



ms 



with a velocity precision of 15 ms l , achieved detection limits ap proaching a Jupite r mass 



for planets in circul ar orbits and period s shorter than ~10 yr. I Wittenmyer et al.l (120061 ) 
combined the data of I Walker et al.l ( 119951 ) with data from the McDonald Observatory planet 
search to achieve a baseline of more than 20 years for a sample of 31 bright solar-type stars. 
Those authors estimated a 99% detec tion limit of 2.0±1.1 M j„ p for planets in Jupiter-like 
orbits (e = 0.0, a=5.2 AU). Recently, ICumming et al.l ( 120081 ) presented a detailed analysis 
of 8 years of Keck Planet Search data (585 stars with N >10 measurements), resulting in 
a typical detectio n limit of 10 ms " 1 for periods less than the duration of the observations. 



From those data, ICumming et al.l ( 120081 ) estimated a giant planet (0.3-10 M Jup ) frequency 



1 For Jupiter, P =11.86yr and K =12.5 ms l . 
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of 10.5% for orbital periods less than 2000 days. 

The Anglo- Australian Planet Search (AAPS) has been in operation since 1998 January, 
and currently monitors 250 stars. The AAPS has achieved a long-term radial-velocity pre- 
cision of 3 ms" 1 or better since its inception, which is enabling the detection of long-period 
giant planets. These planets have typical velocity amplitudes K<15 ms -1 , and their de- 
tection requires superb long-term velocity stability. Some notable recent AAPS detections 
of such planets include HD 134987c (Jones et al. 2010: P = 13.7yr, M sin i=0.82M Jup , 
K = 9ms- 1 ) and GJ 832b (Bailey et al. 2009: P = 9.4yr, M sin i=0.64M Jup , K = 15ms- 1 ). 
The long-term precision of the AAT/UCLES system has enabled the AAPS to be relatively 
efficient at detecting long-period planets. Of the planets discovered by the AAPS, 37±11% 
have periods longer than 1000 days. For comparison, this figure is 29±5% for the Lick & Keck 
program, 18±4% for HARPS+CORALIE+ELODIE, and 14±6% for all other planet-search 
teams. These figures are obtained by counting the number of planets with P >1000 days 
discovered by each group, then dividing by the total number of planets discovered by that 
groupEI As the AAPS now spans 12 years, it is important to make quantitative estimates of 
the population of detectable Jupiter analogs in our sample and in the Solar neighbourhood. 

Here we define a "Jupiter analog" as a planet with a small eccentricity (e < 0.2) and a 
long period (P > 8 yr). That is, a giant planet which plays a dynamical role similar to that of 
our own Jupiter, with a period long enough to imply in situ formation, and an eccentricity 
low enough to suggest a benign dynamical history. In this paper, we present detection 
limits for Jupiter analogs from the 12-year AAPS database. Section 2 briefly describes the 
observational data, and Section 3 discusses the techniques used to calculate detection limits. 
In Section 4, we present the results, and in Section 5 we discuss the results in the broader 
context of the population of long-period planets. 



2. Observations 



AAPS Doppler measurements are made with the UCLES echelle spectrograph (jDiego et al 



199ll ). An iodine absorption cell provides wavelength calibration from 5000 to 6200 A. The 
spectrograph point-spread function and wavelength calibration are deri ved from the i odine 



absorption lines e mbedded on every pixel of the spectrum by the cell (IValenti et al.l 11995 



Butler et al.lll996l ). The result is a precision Doppler velocity estimate for each epoch, along 
with an internal uncertainty estimate, which includes the effects of photon-counting uncer- 
tainties, residual errors in the spectrograph PSF model, and variation in the underlying 



2 Planet data obtained from the Exoplanet Data Explorer at exoplanets.org 
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spectrum between the iodine-free template, and epoch spectra observed through the iodine 
cell. All velocities are measured relative to the arbitrary zero-point defined by the template 
observation. 

The AAPS target list contains 254 stars, of which 180 have been observed for more than 
8 years. Since the aim of this work is to place meaningful limits on Jupiter analogs, here 
we only consider those stars which have more than 8 years of data (N = 180). We further 
restrict the sample to those stars which have been observed more than 30 times (N = 123). 
This is because the reliability of the FAP calculation is strongly dependent on the number 
of data points. 

Table [1] summarises the data characteristics for these 123 stars. For those stars with 
long-term trends indicating a distant stellar companion, a linear or quadratic fit was removed 
from the data before subjecting them to the detection-limit procedure. For stars known to 
host a substellar companion, we did a fit for and removed that orbit and then performed the 
detection-limit computations on the residuals. 



3. Computational Methods 

A set of simulations such as these is only as meaningful as the input assumptions and 
parameters. In this section, we give a detailed discussion of our choices for these simulations. 



3.1. The Detection Limit Algorithm 



The detection limits were computed using the method of IWittenmyer et al.l ( 120061 ). In 
brief, we add a Keplerian signal to the existing velocity data, then attempt to recover that 
signal using a Lomb-Scargle periodogram. The mass of the simulated planet is increased 
until 99% of the injected signals are recovered with FAP<0.1%. For a given mass (or equiv- 
alently, a given velocity amplitude) at a given orbital period, we use a grid of 30 values 
of periastron passage T and, for eccentric orbits, 18 values of the periastron argument u 
spaced evenly every 20 degrees. This makes a total of 30 possible orbital configurations for 
simulated planets with e = 0.0, and 540 configurations for those with eccentric orbits. We 
simulated planets over 100 orbital periods ranging from 1000 to 5000 days, evenly spaced 
in the logarithm of the period. For a simulated planet with amplitude K to be considered 
detectable at the 99% level, we must recover 99% of configurations (30/30 for circular orbits, 
535/540 for eccentric orbits). A successful recovery occurs when the injected period is the 
highest peak in the periodogram, has a FAP less than 0.1%, and the recovered period is 
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within a specified tolerance of the injected period; Section 3.2 gives additional details on 
the selection of this tolerance. In addition to the 99% recovery level, we performed these 
simulations at re covery levels of 90, 70 , and 50%, to maintain c onsistency with our previous 
simulation work (10 'Toole et al.ll2009d ; IWittenmyer et al.ll2010l ). Hereafter, we refer to each 
set of simulations by its eccentricity and recovery rate: for example, the set where e = 0.1 
at 90% recovery is referred to as "e01r90." The result of these simulations is, for each star, 
a plot of the K amplitude (or planet mass) recovered in (99%, 90%, 70%, 50%) of trials at 
each orbital period between 1000 and 5000 days. 

Figure [TJ which shows the distribution of FAP for all simulated signals which were 
considered successfully recovered in the e00r99 trials, demonstrates our reasoning in choosing 
a cutoff criterion of FAP<0.1%. Summed over all 123 stars, the histogram includes results 
from 362368 simulated signals. Though the cutoff criterion was FAP<0.1% (0.001), it is 
evident that the vast majority of recovered signals achieved FAP values that are far more 
significant. The distribution has two peaks, at 10" 4 (N ~12000) and < 10~ 9 (N = 66784). 
The latter of these is an artifact, simply representing the integration over a long tail of 
FAP< 10~ 9 . If most trials resulted in FAP> 10 -3 , the cutoff we have imposed would bias 
the derived detection limits. We would see such an effect in Figured] as a pile-up at the least 
significant bin (FAP= 10~ 3 ). That the highest peaks in Figure [TJ lie at FAP levels at least 
an order of magnitude more significant than the imposed cutoff indicates that our choice of 
cutoff has not biased the results. 

In order to focus on potential Jupiter analogs, we only consider simulated planets with 
e = 0.0, 0.1, and 0.2, in keeping with the e < 0.2 definition of "Jupiter analog" described 
in Section 1. These three values of e provide sufficient sampling of the relevant eccentricity 
range because, for long-period planets (P >1000 days), the orbital eccentricity is typically 
determined only to an accuracy of 0.02-0.04. Figure [2] shows a histogram of the uncertainty in 
eccentricity (cr e ) for the 74 published planets with P >1000 days. The peak of the distribution 
lies at a e ~ 0.02-0.04, with a median value of 0.05. We also note that published estimates of 
cr e are often underestimated due to correlati ons between the Keplerian orbital parameters and 
the non-Gauss i anity o f their distributions (jFordl 120051 ; lO 'Toole et al.ll2009al ). In particular, 
O 'Toole et al.l (j2009al ) found that the 99% confidence interval for eccentricity can be 10-50 
times larger than the traditional uncertainty estimates derived from the covariance matrix 
in a least-squares Keplerian fit. This is particularly important when the signal of the planet 
(defined in O'Toole et al. 2009a as K/<jk) is smaller than about 3. For these reasons, we 
have chosen a coarse grid of eccentricities for these simulations. A finer interval would vastly 
increase the computing time required without adding meaningful information. 
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3.2. False Positives and False Negatives 



AH previous implementations of this detection-limit algorithm ( iWittenmyer et al.ll2006l . 



20071 . 120091 . |2010| ) have used the criterion that the recovered period be within 5% of the 
input period of the simulated planet. The reason for this seemingly arbitrary criterion is 
that in the Lomb-Scargle periodogram, spurious peaks can arise due to the sampling of the 
data. Most common in radial-velocity data are the 1-year and 1-month aliases, the former 
due to targets becoming unobservable as they pass behind the Sun, and the latter due to 
telescope scheduling constraints which usually restrict planet-search observations to bright 
lunations. Harmonics of the injected periodicities (e.g. P/2,2 x P, etc.) also often produce 
significant periodogram peaks. Imposing the "correct-period" criterion reduces the effect of 
these features and thus minimises incorrect detections (false positives). For simulations such 
as these, involving a vast number of attempted planetary detections and using automated 
detection c riteria, it is cri t ical to understand and to minimise (or ideally, to eliminate) false 



positives. lO' Toole et al.l ( I2009af ) presented a detailed discussion of the problem of false 



positives in simulations of planet-search data. 

A potential pitfall in establishing detection criteria to minimise false positives is that, if 
those criteria are too stringent, false negatives become important. False negatives (i.e. incor- 
rect rejections) are substantially more difficult to quantify, and hence are fiendishly difficult 
to control. In practice, because scientists rightfully consider incorrect detections to be a far 
more serious problem, much more effort is spent in eliminating false positives than false neg- 
atives. In the detection-limit algorithm described above, we seek to eliminate false positives 
by the imposition of the "correct-period" criterion. The FAP criterion alone (FAP<0.001) 
eliminates some incorrect detections, but since the FAP calculation is dependent on the 
number of data points, even alias periodicities achieve high significance for larger data sets 
(N > 50). We performed some tests with AAPS data (which have N > 100 for many targets), 
removing the correct-period criterion and relying on only the FAP criterion to eliminate false 
detections. Figure |3] shows the effect of different settings for the FAP criterion in this test. 
As expected, requiring a more significant periodogram recovery (i.e. a lower FAP cutoff) 
reduces the false-positive rate. However, there are two important consequences: 1) The 
false-positive rate remains at a nontrivial level even for extremely stringent FAP levels, and 
2) When the FAP cutoff is set to such low levels, the detection-limit result is dominated 
by false negatives: nearly all of the injected trial signals are rejected, and no meaningful 
detection limit is obtained. 

These tests demonstrated that a correct-period constraint is necessary to obtain use- 
ful detection limits and to control the false positive rate. The next question is then: How 
close to the input period does the recovered period need to be, to be considered a success- 
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ful detection? Too loose a constraint would degrade the scientific value of the detection 
limits obtained from this method, whereas too stringent a constraint results in the auto- 
mated rejection of signals which a rational human investigator would consider sufficiently 
close (e.g. P=3.05 versus 3.04 days). The core question is: "What would a human do?" 
One way to approach this problem is to determine what other humans have already done. 
That is, what is the fractional uncertainty in orbital period at which authors have decided 
to publish a planet detection? Radial-velocity planet search teams generally prefer to wait 
until a complete orbital cycle has been observed before publishing a planet discovery. One 
expects the uncertainty in the fitted period (crp) to drop as more orbital cycles are observed. 
We performed an exhaustive literature search, locating the original discovery paper for every 
radial- velocity planet, and tabulating the period (P), period uncertainty (crp), and the total 
time-span of the observations (A T). In this search, we included only those planets: 1) which 
were originally detected by the radial- velocity method, and 2) where a refereed paper is avail- 
able indicating the above-mentioned quantities. The reason for excluding planets discovered 
by transits is that the transit is a special circumstance which enables the orbital period to 
be determined to extremely high precision. That would skew the relation between ap/P 
and A T, which we seek to apply to the radial-velocity database of the Anglo-Australian 
Planet Search (which includes no transiting planets). The results of this search are shown in 
Figure HI One readily evident feature is that no planets have been published with less than 
0.7 cycles of data (vertical dashed line in the figure). Of the 290 planets, only 24 (8.3%) 
were published with less than one orbital cycle of data. The outlier with an u nusually large 
fractional error in period (o~p/P=17.2%) is HD 149143 ( IFischer et al.ll2006l ). That planet 
was first published with only 17 observations, which may account for the large uncertainty 
in its 4-day period. 

We wish to derive a relation between the data span A T (equivalently, the number of 
cycles), and the period error op/P; however, the number statistics of the known planets 
remain poor. Hence, we performed some additional simulations to fill out the A T-opjP 
plane. For each of 210 stars in the AAPS database, we created 2000 simulated radial-velocity 
datasets. Each simulated dataset used the observation times and velocity uncertainties of 
the original observed data. The simulated planets had periods ranging from 5 to 5000 days 
(in a grid with a uniform spacing of 10 days), velocity amplitudes between 20-50 ms _1 , and 
eccentricities randomly drawn from the range 0.0 : 0.5. Noise (jitter) drawn from a Gaussian 
distribution with a width cr = 5 ms 4 was also added to each simulated velocity mea- 



surem ent. We then fit each dataset with a Keplerian model using GaussFit (jJefferys et al. 



19871 ) and recorded the fitted period and its uncertainty. Since the maximum period used 



sometimes significantly exceeded the length of the available data for certain of the targets 
in the AAPS database, those fits often produced unphysical results, and were discarded. 
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Consistent with the characteristics of published planets described above, we also discarded 
all fits where less than 0.7 cycles of data were present. Of the 420000 simulations, 357681 
remained after these cuts, and the results are plotted in Figure [5J We can now use these 
results to make an informed choice for the correct-period constraint. Since the injected trial 
period is assumed not to be known a priori, it is logical to choose the value of (Jp/P which 
includes 99% of the planets simulated in this section. This results in opjP = 27.7%, which 
is the value used for the correct-period constraint in all detection-limit simulations here. 

Finally, we show in Figure [6] examples of a rejected and an accepted recovery using the 
methods described in this section. An input data file (HD 209653) was chosen at random 
for this demonstration. The left panel shows the periodogram resulting from an injected 
signal with P =1000 days, K = 1ms -1 , and e =0.0. As expected with such a weak signal, 
there is no compelling evidence for it in the periodogram; the highest peaks occur at 9.7 
and 365 days, the latter an artefact of the window function (lower panels). The FAP of 
the 9.7-day peak was 0.018, much less significant than the cutoff of 0.001. The right panel 
of Figure [6] shows the periodogram resulting from the same dataset with an injected signal 
having P =1000 days, K = 10m s _1 , and e =0.0. The peak at 1000 days is quite obvious, 
and its FAP is 0.0005, resulting in a successful recovery. 



4. Results and Discussion 

Due to the large number of stars considered here, it is most efficient to present the 
detection-limit results in terms of the radial-velocity amplitude K which the simulations 
indicate was detectable. Since the targets in the AAPS long-term program are typically 
solar- type stars, one can then estimate the mass thresholds by assuming a 1M Q star. Figure [7] 
shows the distribution of K resulting from the e = 0.0 simulations, summed over all stars 
(123 stars; 12300 K values). A s note d in previous work using this detection-limit algorithm 



(jWittenmyer et al.ll2006l . |2009| . 120101 ) . some injected trial periods can result in apparently 
undetectable signals, due to poor time sampling or large velocity scatter in the input data. 
Circular-orbit results are shown in Figure [7J e = 0.1 results are shown in Figure [HI and 
e = 0.2 results are shown in Figured 

For a given star, the detection limit in K is generally constant over the entire range of 
trial periods shorter than the duration of observations. Figure (TU] shows the detection limit 
in K for three representative stars. It is useful, then, to compute the mean value of K as 
the metric of the quality of a star's detection limit. For each star, we compute the mean 
K (K) and its scatter o~k, then exclude any values more than 3a k away from that mean 
and recalculate the mean and its uncertainty. Summed over all stars and all trials, there are 
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125,100 periods which are shorter than the duration of observations, and a total of 110 K 
values which were rejected by this 3ax clipping process. The discrepant K values tend to 
occur at longer periods, especially those longer than the observations. The results are given 
in Table [2] for all eccentricities and recovery rates. There is a wide range of K, spanning 
two orders of magnitude. This highlights the importance of considering the detection limits 
on a star-by-star basis. For example, HD 19632, the star for which we obtain ed the worst 
detec tion limits (K > 100m s -1 ), has a chromospheric activity index log R' HK (INoyes et al. 
19841 ) of —4.38. This is one of the most active stars on the AAPS target list; "quiet" 



planet-search target stars typically have log R'h k ~ —5.0. Radial-velocity modula tion due 
to starspots can produce such high levels of jitter (IPaulson et al.ll2004t IWrightll2005f ). making 
the determination of meaningful limits extremely difficult. The velocity scatter for the 33 
data points on HD 19632 is 26.2m s -1 , and the data show a strong periodicity near 6 days, 
consistent with its estimated rotation period of 12 days (G. Henry, personal communication). 

The detection limit achievable from a set of radial-velocity data depends largely on two 
factors: the number of data points, and their RMS velocity scatter about the mean. Using 
these results, we can determine an empirical relation between these quantities and the radial- 
velocity amplitude K which can be detected. Such a relation can then be used to estimate 
the amount of observing time required to obtain a robust detection (or non-detection) of a 
particular class of planet. Figure [11] plots the mean K obtained for each star (e = 0.0, 99% 
recovery) versus the quantity RMS/\fN. The 119 stars for which K < 50 ms -1 were used 
in a linear fit, yielding the following relation: 



K =-0.02 + 12.3 (^Ml) ms -i. (1) 

Using this fit, we can make estimates of the number of additional observations required 
to place robust constraints on planets in the AAPS program. For example, a star with 40 
observations at a total velocity rms of 3 ms -1 would yield a detection limit of 5.8 ms" , 
corresponding to 0.5Mj up (P = 12 yr), and 0.1Mj up (P = 50 d). It is important to note that 
the cadence of the observations can also have an impact on the types of planets which are or 
are not detectable. In particular, high-cadence observations are more effective at detecting 
short-period, low-amplitude planets, as discussed further in the next subsection. We note 
that recent modifications to the AAPS observing strategy, such as 20-minute integrations to 
average over stel lar oscillation noise, have resulted in velocity rms of 1-2 ms" 1 s ince 2005 for 
solar-type stars (lO'Toole et al.ll2009b[ Ivogt et al.llioiol : Iwittenmver et al.lfioioh . The AAPS 
plans to obtain a further 6 years of data applying these strategies, which should markedly 
improve the velocity rms achieved (cf. Tabled]), directly resulting in even tighter constraints 
on Jupiter analogs. 



- 10 - 



4.1. Known planet hosts 



Of the 123 AAPS program targets which (1) have more than 8 years of data, and (2) 
have more than 30 observations, 25 are known to host at least one planet (plus one brown 
dwarf host: HD 164427). The planetary parameters are listed in Table A useful sanity 
check is to ask: "Do the detection-limit results indicate that we could have detected the 
(known) planet"? That is, we first removed the known planet's orbit, then added simulated 
planetary signals to estimate the detectability of a given signal. Applying this method to 
the planet hosts, then, serves as a simple check of the degree to which this detection-limit 
method can be trusted. By comparing the velocity amplitude K of the known planets in 
Table [3] to the K detectable for those stars (Table [2]), we find that all but four of those 
planets pass this test at the 99% recovery level: K < K p i anet . The four exceptions are: 
HD 4308b, HD 16417b, HD 2312 7b, and GJ 832b For H P 4 308 and HD 16417, thi s result is 
easily understood by referring to lO 'Toole et al.l (l2009bl ) and lO' Toole et al.l ( 12009a ). Both of 
these low-amplitude, short-period planets were only detectable in AAT data during 48-night 
continuous observing blocks, rather than the several preceding years of more widely-spaced 
observations. Furthermore, these short-period (P < 20 days) planets are wholly different 
from the Jupiter analogs (P > 3000 days) which are the focus of this work. For GJ 832 
and HD 23127, we show the detection limit expressed in terms of K as a function of orbital 
period in Figure H2J In each panel, the known planet is plotted as a large point with error 
bars. We see that for GJ 832, the detection limit at the specific period of the planet is 99%, 
i.e. the planet has a detectability of 99%, though the mean K indicated a detectability of 
only ~90% (averaged over all periods). H ence, these results are consistent with our robust 
detection of GJ 832b feailev et al.l l2009h . For HD 23127 (right panel of Figure Gd, the 



planet lies on the 70% recovery contour; this exam ple illustrates that t 



are more conservative than a human investigator (jO'Toole et alj 



of our simulations also appl ied to the 3 planets orbiting 61 Vir (jVogt et al.ll2010l ). as detailed 



re automated criteria 



2007). This characteristic 



in 



Wittenmver et all fl2010h . 



4.2. Detectability of Jupiter analogs 

In this work, we have adopted the definition of "Jupiter analog" as a giant planet with 
a long period (P^8 yr). By this definition, the AAPS sample includes 3 published Jupiter 
analogs (HD 134987c, GJ 832b, and HD 16069ld). So, to first order, we have a Jupiter- 



3 In this paper, we adopt the designation "c" for the outermost planet in the HD 160691 system, after 
McCarthy etall (|2004h 
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analog frequency of 3/123 = 2.4%. However, this simple calculation assumes that such 
planets are perfectly detectable for all stars in the sample, which is patently false. As we 
have computed circular-orbit detection limits at five recovery levels (99, 90, 70, 50, and 10%), 
we can use these results to apply a rudimentary cor r ection for the relative detectabilities 
(i.e. the completeness) for each star. ICumming et al.l ( 120081 ) detailed a similar technique to 
address the problem of incompleteness for the Keck planet search data. Here, we adopt a 
simple approach and define the survey completeness for a given radial-velocity amplitude K 
and period P as: 



-i iy stars 

f e (P it Ki) = Yl fRj( p i> K i)> ( 2 ) 

stars ■ -. 
3=1 

where fp(P, K) is the recovery rate as a function of K at period P, and N stars is the total 
number of stars in the sample (N = 123). In this way, we account for the detectabilities for 
each star individually, at each of the 100 trial periods. We use the specific detection limit 
Kp obtained for each period from the simulations described above, thus generating five pairs 
of (Kp, recovery fraction). Then, we generate Jr(P, K) for each star by performing a linear 
interpolation between the five pairs of (Kp, recovery fraction). We can then estimate the 
recovery fraction fp(P, K) for any P and K. As an example, suppose we choose K = 10 ms 4 
and P = 1000 days. For HD 10700, which is extremely stable and well-observed, and has 
very tight detection limits (K e= o, = 3.7ms -1 ), a signal of 10 ms" 1 would always be detected, 
and hence fp(1000days, 10ms -1 ) for HD 10700 is 1.0. For a star with poorer detection limits 
such as HD 109200 (K e=Q . = 12.3ms -1 ), we obtain f R (1000days, 10ms -1 ) = 0.795. We can 
see from these examples that if all stars in the sample were stable and well-observed (i.e. if 
selection effects, observing windows, and velocity jitter were unimportant), then every star 
would contribute 1.0 to the sum in Equation (2), giving a survey completeness of 1.0 (100%). 
We could then obtain the planet frequency simply by dividing the number of detections by 
the total number of stars. However, these effects are extremely important for long-term 
radial-velocity surveys, and so we use Equation (2) to obtain a more realistic estimate of 
the completeness of our entire sample as a function of orbital period. Those results are 
shown in Figure [TBI We emphasize that we have not included unpublished planet candidates 
in our estimate of the frequency of Jupiter analogs. In Figure [131 we have summed over 
101 stars, excluding 22 stars which exhibited an unusual artefact arising from the "correct- 
period" criterion discussed in Section 3.2. Consider a data set in which an injected signal 
with P in ~ 3500 days results in P out = 5000 days being recovered by the periodogram. For a 
recovered signal to be accepted, it must be within 27.7% of the injected periodicity. So, for 
P in = 3612d, P out = 5000d is rejected, but at P in = 3671d, P out = 5000d is accepted. This 
resulted in signals which were "undetectable" even at large K becoming eminently detectable 



once Pi n was within 27.7% of 5000d. 

The recovery fraction fn(P, K) in Equation (2) can be used to derive a completeness 
correction for the published detections of Jupiter analogs in the AAPS sample. For each of 
the three stars which hosts a Jupiter analog, we can compute Jr(P, K) at the specific values 
of P and K for that known planet. As we have computed detection limits for eccentricities 
of 0.0, 0.1, and 0.2, we used the results from the eccentricity closest to that of each detected 
planet (GJ 832b: e = 0.2, HD 160691c: e = 0.0, HD 134987c: e = 0.1). This gives 
the following results: GJ 832 - 0.983; HD 160691 - 1.000; HD 134987 - 1.000. Then, we 
compute the survey completeness f c (P,K) to account for the detectability of these specific 
planets around the 120 remaining stars that do not host a Jupiter analog. Here, f c (P, K) is 
computed at the specific values of P and K for each Jupiter analog; those results are then: 
GJ 832 - 0.685; HD 160691 - 0.863; HD 134987 - 0.440. The frequency of Jupiter analogs 
based on this sample, corrected for completeness (detectability), is then given by 



Here, N stars = 123 total stars in the sample, N hosts = 3 which host a Jupiter analog, and 
fn(Pi, Ki) refers to the recovery fractions listed above. In addition, f c (Pi, Ki) is summed 
over the 120 stars which did not host a Jupiter analog, to account for how detectable the 
three found planets would have been around the remaining stars in the sample. 

We can also use the non-detections to compute an upper bound on the frequency of 
Jupiter analogs in the AAPS sample. Using the recovery fraction Jr(P, K) determined as 
above for each star at every trial period, we compute the mean of fn(P, K) over the period 
range 3000-5000 days. Thus, each of the 120 stars which does not host a Jupiter analog has a 
mean recovery fraction fn, with an uncertainty equal to the standard deviation in fn(P, K) 
about that mean. An upper bound on the frequency of Jupiter analogs can then be given by 



where N nonhosts is the number of stars without a Jupiter analog. The result is an upper bound 
on the frequency of such planets with K greater than the selected value. Those results are 
given in Table HI 

Our upper bound of 37% for planets with K > 10 ms 4 and e = 0.0 in the period 
range 3000-5000 days (3 AU< a <6 AU) is consistent with the core-accretion simulations of 




(3) 
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Liu et al.l ( 120091 ). in which 83/220 (37%) of systems resulted in at least one planet matching 



these characteristics. However, our observed frequency of such planets in this sampl e, 3.3%, is 



signifi cantly smaller. This difference arises from the simplifying assumptions in the lLiu et al. 



( 120091 ) simulations. First, each simulated system started with a disk mass of 1 minimum- 
mass solar nebula. A distribution of disk masses may better approximate the formation 
environments of real planetary systems; Liu et al. are performing further tests with this 
modification (Huigen Liu, priv. comm.). Second, the simulations proceeded for 10 7 yr; 
i t is possible that su bsequent dynamical evolution such as planet-planet scattering events 
(IRasio &: Fordl Il996l ) may eject planets in real systems, which would reduce the observed 
occurence rate. 



5. Summary and Conclusions 



Lineweaver fc Gretherl (120031 ) addressed the frequency of Jupiter-like planets, defining 
"Jupiter-like" as those planets with M Saturn < Msin i < 3M Jup and 3 < a < 9 AU. Based on 
the results of eight different Doppler planet surveys, they estimated that 5±2% of Sun-like 
stars host Jupiter-like pl anets. Our estimate of 3 .3% is consistent with theirs. We note that 
this work, unlike that of iLineweaver fc Gretherl (120031 ). is based on detailed simulations on 
a star-by-star basis. This approach, while computationally intensive, is critic al for a proper 



characterisation of the selectio n effects present in radial-velocity data sets (jO'Toole et al. 
2009cl ; Iwittenmver et allboioh . 



The comprehensive study on detectabilities from the Keck planet search given in lCumming et al 



(120081 ) focused on planets with periods less than 2000 days (a ~3 AU for a 1 M Q star). Those 
authors then extrapolated the simulations results to a=20 AU (P=89 years). From their 
Table 2, which gives the giant planet (M p >0.3 Mj up ) occurrence rates for flat and power-law 
extrapolations, they estimated that 2.7±0.8% of stars host planets between 3-6 AU. In that 
range, our m easured occurrence rat e is 3.3±1.4%. Our results are thus entirely consistent 



with those of ICumming et al.l (120081 ) 



In conclusion, we have performed detailed star-by-star simulations on AAPS data with 
a time coverage of 12 years, including the effects of eccentricity, in order to make a robust 
estimate of the frequency of Jupiter analogs in the Solar neighbourhood. Based on our AAPS 
sample, we calculate that no less than 3.3% and no more than 37.2% of stars host a gas giant 
planet in a circular orbit between 3-6 AU. We have also performed a comprehensive analysis 
of the automated detection criteria employed in our simulation method. Salient points arising 
from this work are (1) the uncertainty in the orbital period of a planet increaes dramatically 
when less than one cycle has been observed, (2) no planet has been published with less than 
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0.7 cycles of data, and (3) simulating a large number of planet detections reveals that for 
99% of planets, the fractional uncertainty in period upjP is 27.7%. This figure is useful 
for simulations in which one must determine whether an injected signal has been accurately 
redetected. 
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Table 1. Summary of Radial- Velocity Data 



Star 


N 


RMS 






(ms- 1 ) 


HD 142 


65 


17.6 3 


HD 1581 


87 


3.5 


HD 2039 


42 


15. 7 3 


HD 2151 


172 


4.2 


HD 3823 


64 


5.9 


HD 4308 


101 


4.3 3 


HD 7570 


36 


6.5 


HD 10180 


31 


6.7 


HD 10360 


58 


4.6 1 


HD 10361 


57 


4.3 1 


HD 10647 


35 


14.1 


HD 10700 


209 


3.6 


HD 11112 


30 


8.3 1 


HD 13445 


46 


21. 1 1 ' 3 


HD 16417 


103 


3.6 3 


HD 19632 


33 


26.2 


HD 20766 


37 


10.9 


HD 20782 


40 


5.8 3 


HD 20794 


121 


3.2 


HD 20807 


81 


4.3 


HD 23127 


38 


12. 6 3 


HD 27442 


74 


7.3 3 


HD 28255A 


59 


7.3 1 


HD 28255B 


39 


24.1 1 


HD 38382 


32 


5.1 


HD 38973 


30 


4.7 


HD 39091 


54 


5.6 3 


HD 43834 


115 


5.1 


HD 44594 


30 


5.8 


HD 45701 


30 


5.9 2 


HD 53705 


115 


4.5 


HD 53706 


32 


2.9 


HD 55693 


31 


7.2 


HD 59468 


33 


5.2 


HD 65907a 


53 


6.5 


HD 70642 


37 


4.8 3 


HD 70889 


34 


17.6 


HD 73121 


34 


6.2 


HD 73524 


76 


5.2 


HD 73526 


33 


8.0 3 


HD 74868 


33 


7.8 


HD 75289 


38 


5.9 3 


HD 76700 


40 


6.4 3 


HD 78429 


32 


8.9 


HD 84117 


113 


5.5 
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Table 1 — Continued 



Star 


N 


RMS 






(ms- 1 ) 



HD 88742 


30 


12.8 


HD 92987 


44 


5.9 2 


HD 93385 


36 


6.1 


HD 96423 


33 


5.7 


HD 101959 


36 


6.4 


HD 102117 


53 


4.4 3 


HD 102365 


137 


3.2 


HD 102438 


42 


4.3 


HD 105328 


36 


5.9 


HD 106453 


36 


11.2 


HD 107692 


35 


11.6 


HD 108147 


55 


19.0 


HD 108309 


49 


3.5 


HD 109200 


30 


4.6 


HD 114613 


188 


5.7 


HD 114853 


43 


5.5 


HD 117618 


65 


5.2 3 


HD 120237 


42 


10.5 


HD 122862 


88 


4.4 


HD 125072 


64 


5.0 


HD 128620 


89 


3.5 2 


HD 128621 


126 


3.4 2 


HD 129060 


37 


37.6 


HD 134060 


82 


5.5 


HD 134330 


34 


5.7 


HD 134331 


48 


5.3 1 


HD 134606 


46 


5.3 


HD 134987 


60 


2.5 3 


HD 136352 


134 


4.6 


HD 140901 


76 


10.3 


HD 144628 


42 


4.0 


HD 147722 


57 


16.8 


HD 147723 


60 


9.1 


HD 154577 


31 


4.4 


HD 155974 


39 


7.5 


HD 156274b 


88 


6.4 1 


HD 160691 


155 


2.4 3 


HD 161612 


37 


3.7 


HD 164427 


42 


6.0 3 


HD 168871 


56 


4.9 


HD 177565 


82 


4.0 


HD 179949 


61 


11. 1 3 


HD 181428 


34 


7.9 


HD 183877 


32 


5.3 


HD 187085 


49 


5.7 3 
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Table 1 — Continued 



Star 


N 


RMS 






(ms- 1 ) 



HD 


189567 


73 


5.8 


HD 


190248 


187 


4.1 


HD 


191408 


153 


4.2 1 


HD 


192310 


133 


4.0 


HD 


192865 


37 


11.0 


HD 


193193 


44 


5.8 


HD 


193307 


71 


4.2 


HD 


194640 


64 


4.8 


HD 


196050 


49 


7.8 3 


HD 


196378 


33 


6.7 


HD 


199190 


41 


3.8 


HD 


199288 


56 


5.3 


HD 


199509 


30 


6.6 1 


HD 


202560 


34 


4.9 


HD 


204385 


34 


6.8 


HD 


204961 


32 


5.7 3 


HD 


207129 


100 


4.9 


HD 


208487 


41 


6.0 3 


HD 


208998 


32 


8.4 


HD 


209653 


30 


4.4 


HD 


210918 


57 


4.6 


HD 


211317 


38 


4.3 


HD 


212168 


34 


4.9 


HD 


212708 


34 


4.1 1 


HD 


213240 


33 


4.5 3 


HD 


214953 


72 


6.8 


HD 


216435 


69 


7.0 3 


HD 


216437 


46 


4.9 3 


HD 


217958 


30 


9.4 1 


HD 


217987 


33 


8.9 


HD 


219077 


57 


5.4 2 


HD 


221420 


66 


4.8 1 


HD 


223171 


54 


6.3 



1 Residuals to a linear fit. 

2 Residuals to a quadratic fit. 

3 Residuals after removal of 
known planet (s) orbit. 
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Table 2. Summary of Detection Limits 



Star 


Recovery Rate 
(percent) 


Mean K Detectable (m i 
e = 0.0 e = 0.1 


e = 0.2 


HD 142 


99 


24.8±2.4 


25.9±2.8 


27.8±4.6 




90 


23. Oil. 5 


24.1±2.1 


25.1±2.1 




70 


18.8±2.4 


20.5±2.5 


21.7±10.3 




50 


15.0±2.2 


16.4±1.5 


16.7±2.5 


HD 1581 


99 


7.2±9.0 


6.4±5.4 


6.0±3.3 




90 


4.5±1.4 


4.9±2.1 


5.7±3.6 




70 


3.1±0.6 


3.5±0.5 


3.6±0.8 




50 


2.6±1.1 


2.8±1.1 


3.0±1.4 


HD 2039 


99 


42.6±11.5 


45.9±13.1 


51.0±13.7 




90 


39.2±11.0 


39.5±9.7 


43.8±15.2 




70 


29.9±13.7 


34.6±24.4 


31.0±9.7 




50 


25.0±9.1 


24.4±7.4 


27.0±9.6 


HD 2151 


99 


6.9±2.9 


7.0±1.6 


8.0±2.6 




90 


6. Oil. 2 


6.2±1.1 


6.4±0.8 




70 


5.2±0.8 


5.3±0.6 


5.5±0.6 




50 


4.9±1.4 


4.7±0.7 


4.8±0.7 


HD 3823 


99 


10.4±3.3 


12.2±6.5 


12.8±4.9 




90 


9.8±4.6 


10.2±3.4 


10.9±4.7 




70 


6.5±1.8 


7.3±0.9 


7.7±0.9 




50 


4.8±2.2 


5.0±2.2 


5.2±2.3 


HD 4308 


99 


9.1±8.3 


9.4±7.7 


10.8±8.9 




90 


9.4±9.3 


8.0±5.2 


9.0±6.8 




70 


7.0±3.8 


7.3±4.2 


7.1±3.7 




50 


4.9±2.1 


5.2±2.1 


5.2±1.9 


HD 7570 


99 


13.3±1.2 


14.2±1.9 


15.6±3.5 




90 


11.6±0.8 


12.6±1.1 


13.1±1.2 




70 


9.7±1.0 


10.4±0.5 


10.9±0.7 




50 


8.3±0.9 


8.7±0.7 


9.1±0.8 


HD 10180 


99 


16.3±3.1 


17.2±2.3 


19.8±3.4 




90 


14.7±2.3 


15.5±2.7 


17.2±3.6 




70 


12.1±1.7 


12.8±1.7 


14.0±2.7 




50 


9.9±2.0 


11.0±1.6 


11.3±2.2 


HD 10360 


99 


9.5±5.0 


9.4±3.9 


10.0±3.8 




90 


7.7±1.0 


8.3±2.3 


9.4±6.5 




70 


5.9±0.6 


6.3±0.8 


6.5±0.8 




50 


4.5±0.7 


4.8±0.6 


5.0±0.7 


HD 10361 


99 


7.2±1.2 


7.5±0.8 


7.8±0.9 




90 


6.6±0.5 


6.8±0.6 


7.1±0.7 




70 


5.6±0.5 


5.8±0.5 


6.0±0.5 




50 


4.9±0.5 


5.8±0.5 


5.2±0.4 


HD 10647 


99 


37.3±7.9 


41.6±10.9 


43.3±11.2 




90 


36.0±9.0 


37.6±9.9 


40.5±14.7 




70 


28.3±4.2 


29.9±6.1 


32.5±7.2 




50 


23.6±4.5 


24.5±4.8 


26.0±5.2 


HD 10700 


99 


3.7±0.4 


4.6±2.9 


4.2±0.5 
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Star 


Recovery Rate 
(percent) 


Mean K Detectable (hie 
e = 0.0 e = 0.1 


- 1 ) 

e = 0.2 




90 


3.4i0.4 


3.6i0.4 


3.8i0.5 




70 


o 1 _l_n q 
z. l±U.o 


1 _l_n q 
z. l±U.o 


2.3i0.8 




50 


l.OztU.O 


1.0±U.O 


1.6i0.5 


HD 11112 


99 


OO 74-0 (\ 
ZZ. / ItZ.O 


OA y-i-Q 1 

ZD. / ±y.i 


28.4±2.9 




90 


on /i_l_i n 
zU. 4±l.y 


OO 1 ^0 c; 
ZZ. l±Z.O 


25.4±4.6 




70 


1 1 e:_l_o n 
1 / .oztz.U 


lo.0±l.D 


20.0±2.3 




50 


1 A Q-l-O A 

14.o±z.4 


i c; "7-1-1 1 
lo. / ±z.l 


16.8±2.3 


T TT~\ 1 O A t f 

HD 13445 


99 


A A A _1_1 O A 

41.4±lz.4 


44.z±lo.4 


48.2±14.0 




90 


Ol .Z± / .0 


07 o^fi 1 


41.8±8.1 




70 


QQ 1 

oo. Izto.o 


QO 1 ^/1 O 

oz. 1±4.Z 


36.5±7.5 




50 


O/l 1 i 
z4. 1±D. 1 


OO c^^fi 1 
ZZ. 0±0.1 


23.7±6.3 


T TT"\ -I f A ~t r7 

HD 16417 


99 


O.Uztl.o 


K £4-0 O 
0.0±Z.Z 


6.1±2.1 




90 


4. / ±1.0 


A QJ-A A 

4.o±1.4 


5. Oil. 4 




70 


o. t ±u.y 


o.o±U.y 


3.9±0.8 




50 


o /i _l_n n 


Q_l-i n 


2.8±0.9 


HD 19632 


99 


lz4.1±oo.o 


1 Ofi O-LQfi O 
lZD.Z±OD.Z 


129.8±36.3 




90 


lUo.oitzo.o 


111. / ±00. Z 


123.3±48.4 




70 


OD.Oltz/ .0 


oo.o±oU.D 


88.7±23.3 




50 


70 n_l_oo a 
I o.U±zz.4 


I .U±o0.4 


86.9±29.7 


TTT"\ *w i T / ' / ' 

HD 20766 


99 


00 /l_l_o 1 
Zo.4ztz. 1 


O/l O^O /I 
Z4.y±Z.4 


27.6±3.0 




90 


Ol O-l-O A 

zl.z±z.4 


00 /t_j-i n 
zz. 4±i.y 


24.0±2.2 




70 


1 O Q_l_Q K 

Iz.ozto.o 


1 O *7 

Iz. ( ±0. / 


13.4±9.3 




50 


y.U±D.o 


y.o± / .1 


9.9±7.6 


HD 20782 


99 


1 1 n_l_i n 
ll.Uztl.U 


1 n^o q 


12.7±1.1 




90 


in 1 _l_n 7 
lU.lztU. / 


lU.D±U.y 


11.2±1.2 




70 


0. / ±U.D 


1 ^n 


9.5±0.6 




50 


7.2il.O 


7.8±0.4 


8.3±0.6 


ttt\ r\r\ r ~7/~\ A 

HD 20794 


99 


3.9il.6 


4.2±2.2 


4. Oil. 5 




90 


3.3i0.7 


3.6±0.9 


3.9il.4 




70 


2.6i0.3 


2.7±0.3 


2.8i0.4 




50 


2.3i0.5 


2.3±0.4 


2.6il.O 


HD 20807 


99 


8.9±5.0 


8.7±3.2 


9.2i2.9 




90 


7.0±2.0 


7.6±2.4 


7.8i2.2 




70 


4.8il.l 


4.9±0.9 


5.2il.l 




50 


3.3il.l 


3.3±1.1 


3.3il.2 


HD 23127 


99 


27.4±2.4 


29.4±3.8 


32.2i3.9 




90 


25.5±2.0 


26.3±2.1 


28.4i2.5 




70 


22.8il.8 


23.3±1.7 


24.4il.8 




50 


20.7il.5 


21.2±1.4 


22.2il.4 


HD 27442 


99 


10.0±3.8 


9.9±1.2 


ll.lil.7 




90 


8.7il.O 


9.1±1.2 


9.6il.l 




70 


7.4i0.9 


7.8±0.7 


8.1il.l 




50 


6.2i0.9 


6.4±0.8 


6.7i0.8 


HD 28255A 


99 


17.5±7.1 


18.3±6.8 


21.0il0.6 




90 


14.0±3.0 


15.0i4.2 


36.9i21.0 
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Star Recovery Rate Mean K Detectable (m: 
(percent) e = 0.0 e = 0.1 


e = 0.2 




I u 


11.4±3.5 


11.7i2.0 


11.8±2.0 




CO 

ou 


q n_i_Q p. 


y.4±Z.O 


9.8±2.7 


liU ZoZOotS 


no 

yy 


c;fi A _L1 1 o 


oy.o±iz.o 


65.8±14.7 




on 
yu 


X.A 7 -LI 1 1 
04. i ztlo.l 


£Q c;4.1(1 


58.3±12.0 




f U 


4 I .Olt I .O 


Kr\ O-L 11/1 


53.8±16.7 




KO 
OU 


A O 1 _L7 


/I *7 1 Xt^ Q 

4 / . l±0.o 


47.0±8.0 


un qoqoo 
rlU ooooZ 


on 

yy 


1 O A _LQ Q 

lU.4±o.o 


1 O O-Lft £ 


12.8±4.9 




yu 


O QJ-A a 


1 n O-LQ /I 


10.9±4.7 




i u 


D.Oltl.O 


V Q-l-fl O 

( .oztu.y 


7.7±0.9 




ko 
OU 


/I Q-LO O 

4.o±Z.z 


k n_l_o o 


5.2±2.3 


TJT\ QQ07Q 

liJJ ooy / o 


no 

yy 


lo.o±4. 1 


i c: a J-O o 

1D.4±Z.Z 


18.2±3.6 




oo 
yu 


ll.Oltl.l 


1 9 ft^1 ft 


14. Oil. 7 




TO 

f u 


O Q-LO 

y.oztu.o 


i n 1 n 


11. Oil. 1 




c;o 
OU 


Q O-LO ft 


o 7_l_n ft 


9.1i0.6 


riD oyuyi 


oo 

yy 


O 7_LQ 7 

y. t ±o. ( 


in k 


10.4i2.1 




oo 

yu 


Q Q_L 1 Q 


Q 7-1-1 ft 
O. / ±1.0 


9.3il.8 




70 
I U 


ft 7-LO O 




7.3il.0 




ou 


4.y±i.0 


K ft_l_1 1 
O.Oztl. 1 


5.9il.2 


liD 4ooo4 


oo 

yy 


o.o±4.0 


n Q_l_ft n 

y.ozto.u 


9.9i6.6 




oo 
yu 


7 1 Q 


/ .y±z.o 


8.5i3.1 




f U 


ft oj-1 c; 


D.Oztl.D 


6.9il.7 




KO 

ou 


A OJ-1 O 
'L.Z±1.Z 


4.0±1.U 


4.8il.O 


liU 44oy4 


oo 

yy 


1 ft OJ-C: c; 


i ft QJ- 1 n 

io.o±i.y 


18.3i2.1 




oo 

yu 


lo. / ±Z.Z 


1 /I C^O ft 


15.4il.5 




f u 


1 1 7-1-1 O 

11. i ztl.U 


lO Oil ft 


13.7i3.0 




ou 


1 O A _LO ft 

1U.4±U.0 


iu.y±u.o 


11.6il.O 


T TTv J r TA "1 

HD 45701 


99 


,1,1 *— > i /\ j \ 
22. 7i9. 9 


24.6±7.1 


26.2i7.5 




90 


19.5±6.1 


22.6±11.1 


21.1i7.0 




70 


15.3±4.2 


16.3±4.1 


17.0i4.4 




50 


12.5±1.9 


13.0±1.9 


13.8i2.1 


HD 53705 


99 


6.0±1.4 


6.3±1.2 


6.8il.l 




90 


5.9±1.7 


6.1±2.3 


6.4il.9 




70 


5.3±1.0 


5.5±1.1 


5.9il.4 




50 


4.4±0.5 


4.6±0.6 


4.8i0.6 


HD 53706 


99 


6.0±0.5 


6.7±0.7 


7.7il.5 




90 


5.5±0.2 


5.7±0.3 


6.2i0.4 




70 


4.9±0.3 


5.1±0.3 


5.4i0.3 




50 


4.4±0.4 


4.6±0.3 


4.8i0.3 


HD 55693 


99 


18.5±2.7 


20.6±3.5 


24.7i6.2 




90 


16.6±2.3 


17.8±2.4 


19.8i3.0 




70 


13.9±1.8 


14.8±1.5 


15.9il.6 




50 


11.8±2.1 


12.5±1.9 


13.3i2.2 


HD 59468 


99 


11.4±1.2 


12.1±1.3 


13.4il.4 




90 


10.1±0.7 


10.6±0.9 


11.4i0.9 




70 


8.4±0.8 


9.1±0.5 


9.7i0.5 
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Star Recovery Rate Mean K Detectable (ms 1 ) 





(percent) 


e — u.u 


o — n i 
e — u.± 


e — u.z 




50 


7.0±1.2 


7.4±1.0 


7.8±1.1 


HD 65907a 


99 


10.1±1.5 


10.8±1.6 


11.8±2.2 




90 


8.8±1.0 


9.6±1.5 


10.1±1.5 




70 


6.0±1.8 


7.1±0.8 


7.5±0.7 




50 


4.3±2.0 


4.4±2.1 


4.6±2.1 


HD 70642 


99 


lO.Oil.O 


11.2±2.0 


12. Oil. 7 




90 


8.8±0.7 


9.5±0.9 


10.3±1.0 




70 


7.4±0.6 


7.8±0.5 


8.2±0.7 




50 


6.5±0.6 


6.8±0.4 


7.2±0.3 


HD 70889 


99 


58.5±25.7 


59.7±21.6 


74.2±29.8 




90 


51.1±17.4 


54.2±21.2 


58.4±20.9 




70 


39.6±7.6 


41.6±9.1 


45.1±10.9 




50 


31.9±8.0 


32.4±6.3 


33.4±5.7 


HD 73121 


99 


12.3±1.1 


13.3±0.9 


15.1±2.7 




90 


11.1±0.5 


11.7±0.7 


12.7±0.7 




70 


9.7±1.0 


10.0±0.6 


10.7±0.8 




50 


8.6±1.2 


8.9±1.0 


9.4±1.0 


HD 73524 


99 


10.7±4.0 


11.7±4.5 


12.2±3.3 




90 


8.8±1.8 


9.3±2.0 


10.3±3.7 




70 


6.9±0.9 


7.3±0.7 


7.8±0.9 




50 


5.2±0.8 


5.5±0.8 


5.7±0.8 


HD 73526 


99 


22.8±6.2 


25.3±8.3 


29.1±11.0 




90 


19.5±3.4 


21.4±5.1 


23.6±6.7 




70 


17.8±3.1 


18.6±4.1 


19.2±3.0 




50 


15.2±1.4 


15.8±1.8 


16.9±2.6 


HD 74868 


99 


18.0±1.6 


19.4±2.1 


21.8±3.1 




90 


16.2±1.3 


17.2±1.4 


18.7±1.7 




70 


12.9±2.2 


14.3±1.4 


14.8±1.7 




50 


9.5±3.6 


10.2±3.4 


10.9±3.6 


HD 75289 


99 


11.4±1.0 


12.3±1.4 


13.8±1.6 




90 


10.5±0.7 


11.1±0.8 


11.9±1.0 




70 


8.7±1.0 


9.6±1.0 


9.8±1.1 




50 


7.3±1.3 


8.0±1.4 


8.4±1.5 


HD 76700 


99 


14.3±2.0 


16.3±4.3 


17.6±3.7 




90 


13.4±1.8 


13.9±2.0 


15.1±2.4 




70 


12.2±4.4 


11.4±2.1 


12.4±2.3 




50 


9.9±3.2 


10.2±3.4 


10.6±3.6 


HD 78429 


99 


13.3±1.2 


14.2±1.9 


15.6±3.5 




90 


11.6±0.8 


12.6±1.1 


13.1±1.2 




70 


9.7±1.0 


10.4±0.5 


10.9±0.7 




50 


8.3±0.9 


8.7±0.7 


9.1±0.8 


HD 84117 


99 


8.4±2.9 


8.9±3.5 


9.1±3.3 




90 


7.4±1.2 


7.8±1.4 


8.3±1.9 




70 


6.4±1.0 


6.8±1.0 


7.2±1.1 




50 


5.8±0.8 


6.0±0.8 


6.3±0.9 



-26 - 



Table 2 — Continued 



Star 


Recovery Rate 
(percent) 


Mean K Detectable (m 
e = 0.0 e = 0.1 


s- 1 ) 
e = 0.2 


HD 88742 


99 


31. Oil. 5 


on Q_l_o ft 


40.5i9.1 




90 


ZO.O±l.l 


on o_l_i er 


32.9i2.6 




70 


ZO.O±O.U 


O/t e;4_o c; 
z4.o±z.0 


26.4±2.9 




50 


1 ft 04- £ 7 


1 7 74- 1 


18.8i5.5 


HD 92987 


99 


O ft4_fl ft 


1 n o4_i n 
lU.zztl.U 


10.5i0.7 




90 


o.y±U.O 


n 1 4_n /l 
y. 1±U.4 


9.7i0.8 




70 


Q 1 -LA A 
0. 1±U.4 


o o_l_n q 


8.8i0.3 




50 


7 o4_n Q 


7 Q_l_n Q 


7.9i0.5 


HD 93385 


99 


1 Q A 4-1 Q 


i4.y±i.o 


17.1i3.0 




90 


1 1 0-1-1 1 

n.y±i.i 


1 O 7-1-1 1 
1Z. ( ±1.1 


13.9il.5 




70 


y.y±l.U 


1 n c;_l_n ft 


11.2il.l 




50 


74-1 K 


Q 1 4-1 Q 
y. l±l.o 


9.5il.3 


HD 96423 


99 


1 o Q_l_i n 


1 Q ir _i_ -1 o 

lo.o±l.z 


16.3i4.7 




90 


1 i ft4_n q 
ll.OztU.o 


in < _i_n e 
lz.4±U.o 


13.2i0.9 




70 


1 n /i _l_n c; 


1 n o4_n a 

iu.y±u.4 


11.5i0.5 




50 


y.o±u.o 


n o4_n a 

y.y±u.4 


10.4i0.4 


HD 101959 


99 


1 Q 04-1 O 


1 A Q4-1 K 
14.O±1.0 


16. Oil. 8 




90 


1 1 o_l_n o 


1 1 n4_n n 

n.y±u.y 


13.7il.2 




70 


lU.O±U. / 


1 a e;4_n 7 
lU.O±U. ( 


10.9i0.6 




50 


o i 4-n ft 
y. 1±U.D 


o e;4-n 7 


9.9i0.6 


HD 102117 


99 


7 n4_n 7 
/ .U±U. / 


7 ft4_n n 


8.5i0.6 




90 


ft ft4_n c; 


ft 74_n k 


7.3i0.6 




70 


ft i 4_n 7 


ft o4_i n 


6.4i0.9 




50 


K A 4_n 7 


/i 4_n 


5.7i0.6 


HD 102365 


99 


q fi4-n 7 


/i 1 4-1 n 

4. 1±1.U 


4.9il.9 




90 


q c;4_n 7 


Q 74_n 7 
O. / ±U. ( 


4.1il.2 




70 


3.1i0.6 


3.4i0.7 


3.5i0.7 




50 


2.4±0.4 


2.6±0.2 


2.7i0.3 


HD 102438 


99 


7.5±0.8 


8. Oil. 3 


9.3i3.9 




90 


6.8±0.8 


7.3±0.8 


7.5il.3 




70 


5.1±0.6 


5.8±0.6 


5.8i0.7 




50 


3.8±0.7 


4.2±0.5 


4.3i0.6 


HD 105328 


99 


11.5±0.9 


12.3il.O 


13.6i2.1 




90 


10.5±0.6 


10.9±0.6 


11.8i0.8 




70 


8.9±0.8 


9.5±0.3 


9.8i0.6 




50 


7.5±1.3 


7.9il.l 


8.2il.2 


HD 106453 


99 


24.1±6.0 


25.5±3.1 


29.0i3.7 




90 


22.1±3.9 


23.2±4.1 


23.8il.4 




70 


19.0±3.0 


19.4±2.1 


21.4i4.4 




50 


17.0±2.7 


18.9±5.6 


20.8i8.7 


HD 107692 


99 


29.2±8.3 


29.0±2.4 


32.5i3.7 




90 


24.8±4.2 


26.8±4.6 


29.5i6.7 




70 


17.8±3.8 


19.7±4.1 


20.5i4.6 




50 


14.8±3.2 


15.4±3.1 


16.4i3.3 


HD 108147 


99 


42.5±10.4 


44.0±8.7 


49.5il0.4 
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Star 


Recovery Rate 
(percent) 


Mean K Detectable (ms 
e = 0.0 e = 0.1 


e = 0.2 




90 


A 1 7_l_ 1 c n 

41. 7 ±15. U 


4U.2±8.1 


42.1±6.5 




70 


A 1 7X00 A 
41. / ±.Zz.4 


o4.0±0.0 


36.0±6.3 




50 


ol.O±D.4 


QQ A J.7 A 


35.6±7.9 


HD 108309 


99 


^ .1 _i_n q 


£ fl-l-1 Q 
u.Oztl.o 


5.7±0.7 




90 


a o i n /i 


k o_l_n k 


5.9±3.1 




70 


/i i i n q 


/t Q_l_n q 
4.o±U.o 


4.3±0.3 




50 


q ^ i n /i 
o.O±U.4 


q 7_l_n /i 
o. / ±U.4 


3.8±0.5 


HD 109200 


99 


lz.o±l.U 


io.Oztl.O 


16.4±4.8 




90 


i n 7J.fi 7 
1U. / ±U. / 


11/1 _l_n 7 
11.4±U. ( 


12.9±1.1 




70 


n ^ i o r 
y.oztU.D 


n Q_l_n /i 


10.6±0.4 




50 


o.oztU.O 


o.0±U.4 


9.3±0.4 


HD 114613 


99 


/ .ozto. / 


7 Q4-/1 1 


8.5±4.3 




90 


ft QJ-Q Q 


7 1 c; 


7.2±3.4 




70 


K A 1 O /I 


O.Oztz.O 


5.6±2.5 




50 


/i /IJ_1 n 

4.4±i.y 


A A -1_ 1 C 


4.6±1.9 


HD 114853 


99 


y.i±u. r 


y.oztU.o 


10.5±0.9 




90 


O.OltU.O 


fi o_i_n k 


9.3±0.3 




70 


/ .OltU.O 


7 o_l_n A 
1 .y±U.4 


8.3±0.5 




50 


ft ft I n c 


7 1 _l_n c; 
( .l±U.O 


7.2±0.7 


HD 117618 


99 


7 Q_|_1 /I 


O.OltZ.O 


10.0±5.5 




90 


7 /IJ_1 /I 
/ .4±1.4 


7 c:_l_n n 

/ .o±u.y 


8.1±1.4 




70 


ft 1 J_1 n 

U.1±1.U 


D.O±l.o 


7.1±1.8 




50 


4. / ±l.o 


4.y±1.4 


5.4±1.4 


HD 120237 


99 


i q n_i_o n 
lo.y±^.U 


on n_l_o /l 


21.6±2.8 




90 


1 7 O-U 1 £ 
1 ( .zztl.u 


1 Q OJ_1 o 

lo.zztl.y 


19.3±2.0 




70 


14.z±l. f 


i c i ii n 


15.7±1.2 




50 


11.9±2.3 


12.5±1.5 


12.8±1.9 


TTT~\ 1 o c\ O / ■ O 

HD 122862 


99 


8.3±6.4 


7.8±3.8 


8.0±3.0 




90 


6.5±2.5 


6.5±1.9 


6.6±1.7 




70 


5.0±0.8 


5.1±0.7 


5.3±0.5 




50 


4.4±0.7 


4.5±0.5 


4.7±0.5 


HD 125072 


99 


8.7±1.2 


9.3±1.3 


9.9±1.3 




90 


7.2±1.9 


8.1±1.0 


8.4±1.0 




70 


5.1±2.1 


5.7±2.0 


5.6±1.9 




50 


3.3±2.0 


3.8±2.0 


3.6±2.0 


HD 128620 


99 


6.4±2.2 


7.5±3.6 


9.2±7.5 




90 


5.2±1.7 


5.4±1.4 


6.0±2.0 




70 


3.6±1.1 


3.9±1.1 


4.0±1.3 




50 


2.3±1.1 


3.1±1.1 


2.7±1.0 


HD 128621 


99 


6.2±1.1 


6.9±1.6 


8.1±3.7 




90 


5.5±0.8 


5.7±1.0 


6.1±1.1 




70 


4.2±0.7 


4.2±0.7 


4.4±0.7 




50 


3.1±1.0 


2.9±1.0 


3.1±0.9 


HD 129060 


99 


98.2±21.5 


104.6±24.3 


116.8±28.4 




90 


97.2±25.9 


93.7±16.2 


101.9±17.3 
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Table 2 — Continued 



Star Recovery Rate Mean K Detectable (ms 1 ) 





(percent) 


& — on 
e — u.u 


e — u.i 


e — u.z 




70 


79.8±14.6 


81.8±16.8 


87.2±17.6 




50 


75.8±24.1 


76.7±21.5 


85.8±26.5 


HD 134060 


99 


9.5±6.3 


9.6±3.5 


10.2±2.9 




90 


7.6±2.6 


7.9±1.7 


8.2±1.3 




70 


5.6±0.8 


6.0±0.7 


6.2±0.7 




50 


4.4±1.1 


4.8±1.0 


4.7±1.1 


HD 134330 


99 


12.9±1.3 


13.9±1.8 


15.2±1.7 




90 


11.6±1.1 


12.3±1.3 


13.3±1.7 




70 


10.1±0.8 


10.5±0.7 


11.0±0.9 




50 


8.3±0.7 


8.6±0.8 


9.4±0.7 


HD 134331 


99 


9.4±4.5 


9.2±0.5 


9.8±0.7 




90 


8.6±2.3 


9.0±2.4 


8.8±0.4 




70 


7.2±0.5 


7.6±0.9 


7.9±1.0 




50 


6.0±0.8 


6.5±0.8 


6.6±0.7 


HD 134606 


99 


8.8±1.8 


9.3±2.3 


9.7±0.8 




90 


7.7±0.7 


8.1±1.0 


8.4±0.5 




70 


6.6±0.7 


7.1±0.6 


7.5±0.8 




50 


5.8±0.6 


6.1±0.6 


6.3±0.6 


HD 134987 


99 


3.8±0.4 


4.1±0.5 


4.4±0.7 




90 


3.4±0.4 


3.6±0.4 


3.8±0.4 




70 


2.8±0.5 


3.0±0.5 


3.1±0.6 




50 


2.3±0.7 


2.4±0.8 


2.6±0.8 


HD 136352 


99 


7.0±3.3 


6.9±1.1 


7.7±1.9 




90 


6.1±1.4 


7.0±3.6 


6.5±0.8 




70 


5.0±0.6 


5.4±0.7 


5.7±0.8 




50 


4.3±0.3 


4.7±0.5 


4.8±0.5 


HD 140901 


99 


19.0±8.2 


18.5±4.6 


20.6±7.1 




90 


17.2±4.5 


18.3±7.4 


18.1±5.2 




70 


13.3±2.3 


13.5±2.4 


14.1±2.6 




50 


9.8±2.8 


9.8±3.6 


10.4±3.1 


HD 144628 


99 


7.1±1.3 


7.3±0.7 


8.1±1.1 




90 


6.5±0.4 


6.7±0.4 


7.1±0.5 




70 


5.6±0.4 


6.0±0.6 


6.1±0.5 




50 


4.7±0.4 


5.0±0.3 


5.2±0.3 


HD 147722 


99 


25.4±2.2 


26.5±2.8 


28.8±3.9 




90 


23.9±1.8 


24.5±2.2 


25.3±2.5 




70 


22.9±10.0 


21.1±1.0 


21.7±1.3 




50 


19.2±2.0 


19.6±1.9 


20.5±2.1 


HD 147723 


99 


15.3±2.1 


15.9±2.4 


17.1±2.6 




90 


14.2±1.9 


14.8±2.1 


15.5±2.5 




70 


11.9±1.2 


12.4±1.2 


13. Oil. 6 




50 


9.9±1.6 


10.5±1.4 


11.1±1.6 


HD 154577 


99 


9.9±0.7 


10.8±1.3 


12.7±3.7 




90 


9.0±0.6 


9.6±0.8 


10.5±1.1 




70 


7.9±0.7 


8.4±0.6 


9.0±0.8 
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Table 2 — Continued 



Star Recovery Rate Mean K Detectable (ms r ) 





(percent) 


& — on 


e — u.± 


— no 

e — u.z 




50 


7.1±0.7 


7.4±0.6 


7.9±0.8 


HD 155974 


99 


14.4±1.0 


15.3±1.2 


16.8±1.7 




90 


13.4±0.8 


13.9±0.8 


14.9±1.1 




70 


12.1±1.3 


12.5±1.2 


13.8±3.5 




50 


10.8±1.2 


11.1±0.9 


11.9±1.8 


HD 156274b 


99 


12.2±9.4 


13.6±10.3 


16.9±17.0 




90 


8.3±2.7 


9.4±3.7 


10.9±6.3 




70 


6.3±1.5 


7.1±1.0 


7.2±1.3 




50 


5.3±1.7 


5.5±1.5 


5.7±1.6 


HD 160691 


99 


2.3±0.3 


2.6±0.9 


2.8±1.1 




90 


2.0±0.2 


2.2±0.2 


2.3±0.4 




70 


1.8±0.4 


1.8±0.2 


1.9±0.2 




50 


1.5±0.2 


1.6±0.2 


1.6±0.4 


HD 161612 


99 


7.1±0.5 


7.5±0.7 


8.6±0.9 




90 


6.7±0.7 


6.9±0.4 


7.2±0.4 




70 


5.8±0.7 


6.4±0.7 


6.5±1.2 




50 


5.3±0.5 


5.5±0.5 


5.7±0.6 


HD 164427 


99 


15.9±8.5 


18.0±17.5 


13.9±3.1 




90 


12.6±3.6 


13.9±4.5 


16.0±6.2 




70 


11.5±6.4 


13.0±6.7 


13.5±8.0 




50 


13.4±10.4 


11.9±5.0 


12.8±6.8 


HD 168871 


99 


10.4±8.4 


10.8±7.5 


11.1±7.1 




90 


8.9±4.4 


9.2±4.9 


8.9±3.4 




70 


5.4±2.1 


6.4±1.7 


6.5±3.1 




50 


3.3±1.6 


3.7±1.7 


3.9±1.9 


HD 177565 


99 


7.9±5.0 


7.1±1.5 


7.6±1.4 




90 


6.6±2.8 


6.5±1.2 


6.9±1.6 




70 


3.6±1.8 


3.7±1.9 


4.1±1.8 




50 


2.3±1.5 


2.4±1.5 


2.5±1.5 


HD 179949 


99 


21.8±11.1 


20.9±6.4 


22.1±4.6 




90 


19.1±5.6 


20.8±9.3 


19.5±2.6 




70 


15.3±2.6 


15.4±2.8 


16.9±2.6 




50 


11.1±2.4 


9.9±2.2 


12.2±2.5 


HD 181428 


99 


18.4±1.8 


19.8±2.1 


22.1±3.1 




90 


16.8±1.3 


17.4±1.2 


18.6±1.3 




70 


14.2±1.1 


15.5±1.8 


15.9±2.5 




50 


12.5±1.1 


12.9±1.0 


13.8±1.3 


HD 183877 


99 


12.5±1.7 


14.0±2.6 


16.2±4.0 




90 


11.6±1.9 


12.1±2.4 


13.6±5.9 




70 


9.9±1.0 


10.6±1.4 


11.4±2.3 




50 


8.6±1.0 


9.1±1.1 


9.7±1.4 


HD 187085 


99 


8.9±0.6 


9.2±0.6 


9.9±1.0 




90 


8.4±0.5 


8.7±0.5 


9.0±0.5 




70 


7.5±0.9 


7.5±0.2 


7.8±0.4 




50 


6.8±0.5 


7.1±0.7 


7.3±0.4 
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Table 2 — Continued 



Star Recovery Rate Mean K Detectable (mi 
(percent) e = 0.0 e = 0.1 


e = 0.2 


rlJJ lyyob/ 


99 


10.5±4.2 


1 n A _l_ Q o 

lU.4±o.z 


11.2±3.1 




nn 

yu 


o.0±z.U 


n n_Lo o 
y.U±.Z.z 


9.5±3.0 




7(1 

f u 


7 n_u i q 


7 ft-U 1 o 


7.8±1.3 




OU 


o.y±i.y 


ft Q 


6.6±1.6 


TTT^v 1 OflO/1 Q 

HD lyUZ4o 


on 

yy 


A n I i /i 

4.y±±.4 


Z. ft I 7 
0.O±.Z. ( 


5.8±2.1 




on 

yu 


/i QJ_n n 
4.o±U.y 


/i c i n n 

4.o±u.y 


4.8±1.2 




f U 


o.O±U.O 


Q Qj_n 7 


4.0±0.6 




c;n 
OU 


o n 1 n 


o..Z±U.O 


3.2±0.4 


liU iyi4Uo 


on 

yy 


a a i n 7 
4.4±U. / 


4. r ±U.o 


5.1±0.8 




nn 

yu 


/i i J_n 
4. l±U.O 


a o_i_n ft 
4..Z±U.0 


4.5±0.6 




f u 


q /i i n q 
o.4±U.o 


Q 7 i n o 
o.0±U.z 


3.7±0.3 




OU 


Z.o±U.o 


z.y±U.z 


3.0±0.3 


TT"Pl 1 OOQ1 n 


on 

yy 


o.4±l.o 


q ft i 1 n 


3.8±2.0 




yu 


Q 1 1 1 ft 

O.l±1.0 


Q OJ_ 1 7 


3.4±1.7 




f u 


O ft 1 1 Q 


Z. i ±l.o 


2.8±1.4 




c;n 

ou 


z.4±l.z 


O /I _L1 Q 


2.5±1.3 


riU lyzoOO 


QO 

yy 


Zo.4±0.y 


Q1 O-Ufi ft 


35.9±11.0 




on 

yu 


Oft QJ_ft O 


07 q_lc ft 


32.5±11.4 




f U 


01 C-LQ A 


OQ n_LQ 1 

Zo.U±o.l 


24.4±3.6 




OU 


I 7 O-LQ n 

I I .y±o.u 


1 O c;-UO 7 

iy .0±Z. / 


"1 n T 1 O A 


liD lyoiyo 


on 

yy 


17 1 _L1 A K 

1 / .1±14.0 


1 Q /111 C C 

lo.4±10.0 


19.5±17.5 




on 

yu 


i cr n_Li n q 
10.U±lU.o 


1 ft 7-L1 /I K 

lb. r ±14.0 


17.4±14.8 




f u 


1 Q 1 1 1 O Q 


lz.o±11.0 


11.0±4.3 




c;n 
OU 


n n_LQ o 


q n_i_o n 


9.9±3.6 


rlU lyooU* 


on 
yy 


n /1 -U 7 ft 

y.4± / .0 


o o-i-A n 
y.z±4.u 


11.0±6.0 




nn 

yu 


7 n i o ft 


ft QJ-1 ft 


6.8±1.4 




70 


5.1±0.8 


5.2±0.8 


5.4±0.9 




50 


4.2±1.0 


4.4±0.8 


4.4±0.9 


HD 194640 


99 


10.4±9.1 


8.9±3.3 


8.7±1.3 




90 


7.4±1.1 


7.5±0.7 


7.9±0.8 




70 


6.5±0.7 


6.6±0.8 


6.9±0.9 




50 


6.2±1.0 


6.3±0.9 


6.5±1.0 


HD 196050 


99 


13.6±1.7 


14.3±2.0 


15.7±1.9 




90 


12.9±1.4 


13.4±1.7 


13.9±1.6 




70 


11.6±1.8 


12.1±1.2 


12.2±1.1 




50 


9.5±1.3 


10.9±1.5 


11.2±3.0 


HD 196378 


99 


15.3±1.1 


16.8±1.6 


18.7±1.1 




90 


14.0±0.9 


15.0±1.1 


16.4±1.7 




70 


12.7±1.9 


13.2±0.8 


14.1±1.1 




50 


12.0±1.7 


12.3±1.6 


13.1±1.9 


HD 199190 


99 


6.8±0.6 


7.2±1.6 


7.7±0.5 




90 


6.4±0.5 


6.6±0.5 


7.1±0.6 




70 


5.9±0.7 


6.1±0.9 


6.1±0.4 




50 


5.5±0.6 


5.7±1.0 


5.7±0.5 


HD 199288 


99 


10.2±1.4 


10.7±1.4 


11.4±1.6 
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Table 2 — Continued 



Star Recovery Rate Mean K Detectable (ms 1 ) 





(percent) 


e — u.u 


o — n i 
e — u.± 


e — u.z 




90 


9.2±1.1 


9.6±1.2 


10.2±1.4 




70 


7.5±1.0 


7.8±0.9 


8.2±1.2 




50 


6.8±1.2 


6.9±1.1 


7.3±1.2 


HD 199509 


99 


25.7±10.6 


29.4±14.7 


33.7±19.8 




90 


23.1±10.0 


22.7±8.7 


23.7±8.8 




70 


15.9±3.3 


16.7±4.2 


21.1±9.3 




50 


13.3±2.4 


13.4±2.3 


15.7±4.7 


HD 202560 


99 


10.6±1.0 


11.4±1.2 


12.4±1.5 




90 


9.7±0.7 


10.0±0.7 


11. Oil. 2 




70 


8.2±0.5 


8.5±0.4 


8.9±0.3 




50 


6.9±0.7 


7.2±0.4 


7.5±0.6 


HD 204385 


99 


14.6±0.7 


15.7±0.8 


17.9±2.2 




90 


13.5±0.7 


14.1±0.5 


15.4±0.8 




70 


11.3±1.8 


11.9±0.9 


12.5±1.2 




50 


9.1±1.7 


9.6±1.6 


10. Oil. 7 


HD 204961 


99 


14.0±4.7 


14.9±1.8 


16.7±2.2 




90 


12.8±1.4 


13.2±1.5 


14.2±1.3 




70 


11.4±0.9 


11.8±0.9 


12.7±1.1 




50 


10.3±0.6 


10.7±0.6 


11.2±0.8 


HD 207129 


99 


8.7±1.2 


9.3±1.4 


10.3±1.7 




90 


7.5±1.1 


8.0±1.1 


8.6±1.1 




70 


5.9±0.9 


6.2±0.8 


6.3±0.7 




50 


4.7±1.0 


4.8±0.5 


5.0±0.7 


HD 208487 


99 


11.1±1.2 


12.1±1.6 


13.9±2.2 




90 


10.0±0.8 


11. Oil. 9 


11.4±1.0 




70 


9.2±0.6 


9.6±0.6 


10.2±0.7 




50 


8.6±0.8 


9.0±1.0 


9.1±0.6 


HD 208998 


99 


18.9±1.9 


20.9±5.0 


24.4±11.3 




90 


17.4±1.2 


18.3±1.6 


20.2±3.0 




70 


16.4±1.6 


17.1±2.5 


18.1±3.0 




50 


15.2±1.3 


15.5±1.7 


16.1±1.4 


HD 209653 


99 


10.6±0.9 


12.7±6.0 


12.9±0.8 




90 


9.8±0.7 


10.2±1.0 


11.4±1.8 




70 


8.8±1.1 


9.1±1.3 


9.7±0.7 




50 


8.3±1.2 


8.5±1.5 


9.0±1.4 


HD 210918 


99 


13.6±11.1 


14.5±12.3 


15.9±9.7 




90 


13.5±11.8 


13.0±9.5 


14.6±11.1 




70 


6.9±3.7 


7.5±4.3 


7.7±4.3 




50 


4.6±3.0 


4.7±3.1 


5.0±3.3 


HD 211317 


99 


10.4±3.7 


10.8±2.6 


11.7±1.8 




90 


8.6±1.4 


9.0±1.1 


9.8±1.1 




70 


7.3±0.8 


7.7±0.7 


8.1±1.2 




50 


6.4±0.9 


6.7±0.8 


7.1±1.0 


HD 212168 


99 


10.6±0.7 


11.7±0.7 


12.9±1.0 




90 


9.7±0.6 


10.2±0.3 


10.9±0.5 
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Star 


Recovery Rate 


Mean K Detectable (ms 


- 1 ) 




(percent) 


e = 0.0 


e = 0.1 


e = 0.2 




70 


8.4±0.5 


9.0±1.0 


9.6±1.4 




50 


7.6±0.5 


7.9±0.6 


8.3±0.8 


HD 212708 


99 


13.1±4.3 


15.6±6.0 


21.9±13.0 




90 


12.4±5.6 


11.5±3.8 


12.6±4.5 




70 


9.6±2.7 


11.7±4.1 


12.7±5.7 




50 


7.6±1.3 


8.1±1.4 


8.6±1.6 


HD 213240 


99 


16.5±7.8 


19.1±9.3 


21.2±9.6 




90 


12.5±3.0 


14.1±4.4 


15.4±4.7 




70 


11.3±3.3 


13.5±5.1 


14.7±6.2 




50 


9.1±1.3 


10.0±2.0 


10.8±2.4 


HD 214953 


99 


19.8±17.1 


18.7±13.1 


18.2±8.6 




90 


11.8±3.1 


14.0±3.7 


13.1±5.4 




70 


9.0±6.5 


8.7±4.2 


8.8±4.3 




50 


7.3±3.3 


7.4±3.4 


7.7±3.5 


HD 216435 


99 


8.5±0.5 


8.9±0.7 


10.1±2.1 




90 


8.3±0.5 


8.4±0.5 


8.7±0.6 




70 


7.7±0.7 


7.7±0.9 


8.2±1.0 




50 


6.9±0.6 


7.0±0.7 


7.3±0.8 


HD 216437 


99 


8.4±0.9 


9. Oil. 3 


10.4±3.1 




90 


7.5±0.6 


8.0±0.6 


8.4±0.8 




70 


6.7±0.8 


7.2±1.8 


7.3±1.0 




50 


6.0±0.6 


6.3±0.5 


6.6±0.6 


HD 217958 


99 


32.3±7.3 


36.0±9.1 


43.8±14.7 




90 


30.7±11.7 


30.0±6.5 


35.3±9.8 




70 


27.6±7.8 


26.8±6.8 


31.9±10.0 




50 


23.3±5.4 


23.4±4.9 


26.2±6.8 


HD 217987 


99 


22.0±3.8 


24.0±4.9 


27.5±6.2 




90 


20.1±2.5 


21.0±3.0 


22.8±3.4 




70 


18.6±3.2 


18.2±1.3 


19.6±1.6 




50 


16.6±1.5 


17.0±1.5 


18.1±2.4 


HD 219077 


99 


13.5±8.2 


14.0±7.5 


15.5±8.8 




90 


9.5±2.8 


10.6±3.4 


11.0±3.9 




70 


6.3±2.2 


7.1±1.7 


7.1±2.4 




50 


4.0±2.4 
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Log False Alarm Probability 

Fig. I. — Distribution of false-alarm probabilities for 362368 successfully recovered signals 
from detection-limit simulations for 123 stars (e = 0.0, 99% recovery). Though the cutoff 
criterion was 0.001, the vast majority of signals were recovered at greater significance, and 
18.4% of signals were recovered with FAP< 10~ 9 . Of all trial signals, 23.0% were rejected 
based on the FAP cutoff criterion. 
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Fig. 2. — Distribution of the published uncertainties in eccentricity for the 74 published long- 
period planets (P >1000 days) based on exoplanet.eu, 2010 September 1. The distribution 
peaks in the range <r e =0.02-0.04; the true unc ertain t y in e ccen tricity may be conside rably 
larger due to non-Gaussianities as described in iFordl (120051 ) and lOToole et al.l (j2009al ). 
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Fig. 3. — Mean false positive rate obtained from tests of the detection-limit routine in which 
the correct-period criterion was removed. Incorrectly recovered periods tend to have higher 
FAPs, and should be excluded, but a nontrivial number of false positives occur even at 
extremely stringent FAP levels. 
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Fig. 4. — Fractional error in orbital period as a function of the number of cycles observed, for 
290 published radial-velocity planets. The vertical dashed line indicates 0.7 cycles of data, 
which appears to be the minimum for publication. 
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Number of Cycles 

Fig. 5. — Fractional error in orbital period as a function of the number of cycles observed, for 
362953 simulated radial- velocity planets. As in Figure HI the vertical dashed line indicates 
0.7 cycles of data. Note that the slope of the relation becomes markedly steeper when less 
than one cycle is available. 
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Fig. 6. — An example of a rejected (left panel) and an accepted (right panel) detection using 
the method and criteria outlined in Section 3. Both periodograms result from the addition 
of a signal at P = 1000 days to the velocity data for HD 209653. The injected signals have 
K = 1 ms- 1 (left) and 10 ms 4 (right). 
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Fig. 7. — Distribution of velocity amplitudes (K) recovered for simulated planets with e = 
0.0, at recovery levels of 99, 90, 70, and 50%. Each panel shows results from all 123 stars 
(12300 trial K values). Only K < SOms^ 1 is shown; this range includes 91.7% of trials. The 
vertical dotted line at K = 12 ms^ 1 indicates the radial-velocity signal of Jupiter. 
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Fig. 8. — Same as Figured but for the e = 0.1 results. Only K < 50 ms 1 is shown; 
range includes 90.4% of trials. 
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Fig. 9. — Same as Figured but for the e = 0.2 results. Only K < 50 ms 1 is shown; 
range includes 88.8% of trials. 
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Fig. 10. — Detection limits computed for three representative stars, at 99% recovery. In 
terms of the velocity amplitude K, the detection limit is essentially constant for a given 
star, for periods shorter than the duration of observations. At periods longer than the 
observational data, the injected test signals cannot be robustly detected by our automated 
criteria, as seen for HD 102117 here (T obs =3961 days). 
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Fig. 11. — Radial- velocity amplitude K that can be detected at the 99% confidence level, 
plotted against RMS/ y/(N), for 119 stars. A linear relation can be fit, with a scatter of 3.6 
ms" 1 . The solid line is a weighted fit, which is dominated by the abundance of points at 
K < 15 ms" 1 with small error bars. 
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Fig. 12.— Detection limits for GJ 832 (left) and HD 23127 (right) at e = 0.2. The known 
planet is plotted as a large point with error bars. At the specific period of each planet, the 
detectability is 99% for GJ 832 and 70% for HD 23127. This shows that the automated 
detection criteria used in this work are more conservative than a human investigator. 
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Fig. 13. — Completeness fraction f c (Pi, Ki) (Equation 2) summed over 101 AAPS stars, as a 
function of orbital period and radial-velocity amplitude K. From bottom to top, the curves 
are for K =10, 20, 30, 40, and 50 ms" 1 . The drop-off in completeness at P ~4300 days 
occurs because that is the maximum duration of observations for any star in the sample. 
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Table 4. Jupiter-Analog Upper Limits from the AAPS Sample 



Velocity Amplitude 
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12.6±1.1 


13.6±1.4 


16.2±1.5 


K > 30 
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18.6±1.5 


K > 20 


18.6±1.1 
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23.8±1.6 


A" > 10 


37.2±1.1 


44.8±1.4 


48.8±1.5 



